Superconductivity in hole-doped Cgo from electronic correlations. 
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We derive a model for the highest occupied molecular orbital band of a Ceo crystal which includes 
on-site electron-electron interactions. The form of the interactions are based on the icosahedral 
symmetry of the Ceo molecule together with a perturbative treatment of an isolated Ceo molecule. 
Using this model we do a mean-field calculation in two dimensions on the [100] surface of the crystal. 
Due to the multi-band nature we find that electron-electron interactions can have a profound effect 
on the density of states as a function of doping. The doping dependence of the transition temperature 
can then be qualitatively different from that expected from simple BCS theory based on the density 
of states from band structure calculations. 
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Superconductivity in Cgo has generally been ascribed 
to a phonon mechanism due to strong electron-phonon 
coupling for some Cgo intramolecular modes [1]. How- 
ever, due to their high energy and the narrow electronic 
band width questions have been raised about the effec- 
tiveness of retardation for reducing the strong Coulomb 
repulsion in these materials [2]. In addition, a number 
of features suggest that these materials are very exotic, 
including Mott insulating behavior and so called "bad 
metal" behavior with resistivities which do not saturate 
at high temperatures [3,4]. Such behavior is reminiscent 
of the high-Tc cuprates where electron correlations are 
generally accepted to play a crucial role. 

Another issue which is not well understood is the vari- 
ation of Tc with doping. The various alkali-doped ma- 
terials have Tc 's that are maximized near the half-filled 
LUMO (lowest unoccupied molecular orbital) band, i.e. 3 
electrons per Cgo, and confined to a narrow doping range 
around this. This variation of Tc with doping does not 
correspond to the density of states (DOS) as given by 
band structure calculations [3]. Again, a clear indicator 
that correlation effects are important. 

Here we study the effects of electronic correlations on 
a crystal of Ceo molecules based on the strong intra- 
molecular electron-electron interactions. Our approach, 
which expands on earlier work by Chakravarty et al. [5] 
(see also [6] for similar and independent ideas), can be 
summarized as follows. We solve the Hubbard model 
on a truncated icosahedron, i.e. a single Cgo molecule, 
in second order perturbation theory in the on-site re- 
pulsion U. We do this for the HOMO (highest occupied 
molecular orbital) states given by diagonalizing the tight- 
binding {U — 0) Hamiltonian. Based on the perturbative 
spectrum we formulate an effective interacting Hamil- 
tonian in terms of holes characterized by orbital angu- 
lar momentum and spin. We then consider a crystal of 
Ceo molecules with nearest neighbor hopping and where 
this effective Hamiltonian for the interactions on a single 
Ceo molecule correspond to on-site interactions. Subse- 
quently we do a standard BCS/Hartree-Fock calculation 
on the [100] surface of an fee crystal using this lattice 
Hamiltonian. The Hartree-Fock calculation on a surface- 



and hole-doped Cgo crystal is a model calculation. Nev- 
ertheless we feel that the method presented as well as the 
qualitative features of the results are relevant also to the 
alkali doped materials. 

We find that Tc is peaked close to three holes and 
strongly suppressed at five holes where the DOS based 
on band structure is maximized. This striking deviation 
from the behavior expected from a BCS calculation based 
on band structure is related to a strong renormalization 
of the DOS due to the interactions. As a signature of the 
strong electron-electron interactions we also find that de- 
pending on the details of the interactions and band struc- 
ture there may be non-magnetic Mott insulating states at 
even integer fillings. Mott insulating behavior is indeed 
seen in alkali doped compounds with a doping of two or 
four electrons per molecule [7,8]. In addition, an equiv- 
alent analysis for the LUMO band gives a pair-binding 
interaction which is roughly 60% of that for the HOMO 
band in the relevant parameter regime, which suggests a 
possibility for higher Tc 's for a hole-doped material. 

Perturbation Theory: Let us start by considering the 
following Hubbard model on a single Ceo molecule, 
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where the only non- vanishing hopping integrals are = 
t for nearest neighbor (nn) hopping on pentagons and 
tij — t' for nn hopping between pentagons. We use 
t = 2eV and t' /t — 1.2 in accordance with the values 
used in [5] and allow U to vary. Values of [/ ~ 5 — 12eV 
have been suggested in the literature [3] . 

We do standard second order perturbation theory in 
Hubbard U. Since the Hamiltonian is spin rotationally 
invariant the states fall into degenerate sets correspond- 
ing to irreducible representations of the icosahedral group 
Ih and spin. The states are well characterized by angu- 
lar momentum and only weakly split by the icosahedral 
symmetry. 

The validity of second order perturbation theory for 
the large-U Hubbard model and the neglect of longer 
range Coulomb interaction for this problem has been un- 
der debate [3]. It has been shown by exact diagonal- 
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ization that for small Hubbard clusters (e.g. the 12- 
site truncated tetrahedron) the second order perturba- 
tion theory is qualitatively correct giving positive pair- 
binding energies for moderately large U [9] . In addition, 
longer ranged repulsions are more effectively screened by 
the metallic environment than the high energy second or- 
der processes that give rise to the attraction [10]. Here, 
we explore the consequences of this model assuming that 
it gives a reasonable estimate of the molecular spectrum 
and the resulting pair attraction. 

Effective Interactions: Using the perturbation theory 
results we derive a set of interactions for the crystal. Al- 
though the perturbative result contains terms to higher 
order in fermion operators, we take the effective Hamil- 
tonian 



H, 



eff = eo-ei 



which acts on a space of one particle statcis in the five- 
dimensional Hu representation of Ih- Here cj.^ creates 
a hole with quantum number k = 1,...,5 and spin a, 
ni„ is tli{^ number operator and alike indices are hence- 
forth summed over. The eo, ei and ti are parameters and 
the Tli^^^p^g are tensors chosen such to make the four- 
fermion terms invariant independently under spin and 
icosahedral symmetry. 

Group theory reveals that there are nine such indepen- 
dent four-fermion terms. These can be derived by con- 
structing all two-fermion terms CkaCip transforming in 
a particular representation of spin and angular momen- 
tum and taking the tracing with their hermitian con- 
jugates. We can write for the product of two fermions 
in the representation H of Ih and spin— 1/2, H ® H = 
{A®G® 2H), + (Ti ® T2 ® G)a, 1/2 O 1/2 = 0^ © 1. 
where s and a mean the symmetric and antisymmetric 
parts of the tensor products and where A,Ti,T2,G and 
H are the 1,3,3,4 and 5 dimensional representations of Ih 
respectively. 

The product of two anticommuting fermion operators 
thus reduces into seven irreducible parts, given by find- 
ing the antisymmetric part of the product of angular 
momentum and spin. We then construct the invariant 
four-fermion operators T% with corresponding coupling 
constants ti, labeled according to from which two-fermion 
operators they are constructed using the composite index 

i = {iA, 0, 0), (if, 2, 0), (G, 4, 0), (H, 4, 0), (Ti, 1, 1), (3) 
(T2,3,1),(G,3,1)}. 

Here {Ri,Li,Si) indicates icosahedral representation, 
corresponding angular momentum in the case of full rota- 
tional symmetry, 0(3), and spin, respectively. The ten- 



sors Ti 



are normalized such that they are pro- 
jection operators into the ith irreducible subspace of the 
two-fermion products. We neglect the additional two in- 
variants that can be constructed from tracing the two 



different realizations of H together since these arc not 
allowed under 0(3) and since the deviance from full ro- 
tational symmetry is small. 

The effective Hamiltonian (2) is then used to match the 
spectrum found from the perturbation theory, eo is given 
by the energy of the neutral molecule (filled HOMO), ei 
by the 10-fold degenerate 1-holc states and the 2-hole 
spectrum is in one-one correspondence to the seven four- 
fermion terms with energies eo — 2ei -|- 2ti, which fixes 
ti =ti{U,U'^) as shown in figure 1. 
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FIG. 1. Coupling constants of tlie effective interactions la- 
beled according to {L,S). The split due to icosahedral sym- 
metry is not explicitly labeled. {U is in units of t = 2eV) 

In spite of its relative simplicity this effective model 

reproduces the trend found in the perturbative calcu- 
lation, namely that for moderately large U states with 
low spin and low "angular momentum" have lower en- 
ergy, i.e. Hund's rule is not valid. To make this state- 
ment more transparent we can consider the conventional 
four-fermion operators n^, S'^ and L^, invariant under 
0(3) X SU{2). To fit these to the normal ordered four- 
fermion operators T* we also need the tensor invariants 
of angular momentum Q^QI (no sum over L) where 
Ql = Ql.ij'^la'^-ja transform as the L =2,3 or 4 rep- 
resentation of 50(3) (a 

fined in (3) and with similar expressions for the other 
operators. Apart from the small split due to icosahedral 
symmetry, which we average over, we get new coupling 
constants as shown in fig. (2). (The L = 2,3 invariants 
are normalized as and due to overcompleteness we 
choose not to include the L = 4 tensor invariant.) 
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FIG. 2. Magnitude of the coupling constants of S^, n^, and 
respectively g, u, A, and the L = 2 and L = 3 tensor in- 
variants. {U is in units of t = 2eV) 

We find that for large U, n^, and S'^ dominate 
the energetics and we will subsequently only use these 
as a minimal model expected to capture the important 
physics of the interactions. 

Lattice Hamiltonian: We can now write down the lat- 
tice Hamiltonian 



(4) 



<fr'> 



ur&, 



where r are the sites of the lattice and < rr' > runs over 
the range of intermolecular hopping. This model is quite 
general and could be used also for electron doped systems 
in two or three dimensions [11]. Here we consider a model 
where the charge is confined to the [100] surface of the 
fee crystal, i.e. a two dimensional system where < rr* > 
correspond to nn on a square lattice. From the pertur- 
bative calculation (fig. 2) we take g,X,u > 0. Note that 
the confinement of electron propagation to the surface 
completely breaks the five-fold degeneracy of the HOMO 
states. This will be manifest in the tight-binding part 
of the Hamiltonian which reflects the symmetry group of 
the surface. 

Typical hopping integrals are of the order of O.leV [12] 
which is comparable to the interactions (fig. 2), the prob- 
lem is in an intermediate coupling regime. Nevertheless, 
we do a standard BCS/Hartree-Fock construction, re- 
placing the Hamiltonian (4) by a reduced non-interacting 
Hamiltonian. We keep only spatially uniform supercon- 
ducting mean-fields bkiap = y < Cf^kaCr,ii3 > and 
mean-fields of the number operators nia = y < 
ctj^^Cf >. Wc can then derive the following effective 
Hamiltonian in momentum (p) space 

HmF= Ho + Hpair + HhF = '^{t{p)kl - I^Skl)cl kcr^,la- 

p 

L,S 

where we included a chemical potential /U. We define the 
components of the order parameter 

^kiaff ~ '^L,S X '^klmna/SfS^mniS , (6) 

i:Li=L,Si=S 

where Vl,s = (| - + (6 - ^^^^)A - u and 

Vl,s > corresponds to attraction. 

Assuming no net magnetization the Hartee-Fock terms 
can be written (no sum over /) 

hHF,ii = --^9ni-XTi+u{ni + 2^nk), (7) 

kjil 



where n; = n;| is the total particle number with an- 
gular momentum component I and r = (3n4 -t- 8715, Ans + 
Hi + n5, 4n2 -1-714-1-715, 3711 + n2 + n3 + n5, Sni + n2 + 
n-3 + 714). (In addition there is one off-diagonal compo- 
nent •\/3A(7i5 — 7i4)(Cp 2aCj),2a +h.c.), which is included in 
our calculations but which will in general be small.) 

The tight-binding part of the Hamiltonian for the sim- 
plest (unidirectional) crystal structure takes the form 



t{p)ki = 
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(8) 



where fp = cos(px) cos.{py) and gp = sm{px) sm{py). 

We have taken hopping parameters from [12]. For 
hole hopping we have, ti = —.107,^2 = .198,^3 = 
.134, i4 = -.032, is = -.170, ti2 = .087,^13 = .073, t23 = 
.162,^45 = .115 eV. The Hamiltonian has the symmetry 
t{p) = ~t{p ~\~ (■"'jO)) implying a symmetric band struc- 
ture around zero energy where all bands will be half-filled 
and there are Van-Hove singularities at zero energy at 
(px,Pj,) = (±7r/2,±7r/2). 

The Hartree-Fock terms (7) have interesting proper- 
ties related to their multiband nature. For positive g 
and u it is energetically favorable to fill up as few bands 
as possible for a given particle number. The term Sp ■ Sp 
gives on-site spin-triplet states with higher energy than 
singlets, so that by putting particles in a single angular 
momentum state the energy can be lowered by exclusion, 
and similarly for the on-site charging energy n^, where 
two particles with the same spin and angular momentum 
cannot occupy the same site. The Lp ■ Lp term on the 
other hand gives rise to an anisotropic attraction between 
the components. 

For positive parameters <?, u and A the L = 0, 5 = 
pairing channel of (5) is strongest and we can expect this 
to dominate. But, since the rotational invariancc of the 
Ceo molecules is broken by the lattice, subdominant or- 
der parameters with non-zero angular momentum appear 
and in general all three 5 = order parameters may be 
non-zero. 

What kind of physics can we expect from this model? 
For large g or u there may be Mott insulating states at 
even integer 27i filling where n bands of angular momen- 
tum states will fill up completely. If u is not very large 
compared to g the insulating state will be non-magnetic 
due to the low energy of molecular singlets. (For the 
regime g, not realized here, there may also be mag- 
netic insulating states at odd integer filling.) The pairing 
terms compete with a putative insulating state due to the 
Hartree-Fock terms so that even for large 5 or A there 
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may be s.c. ground states also at even integer fillings, al- 
though a suppression of Tg is likely due to the low DOS 
when the bands are nearly filled or empty. In general we 
can expect the Hartrcc-Fock terms to completely recast 
the DOS compared to that given by the band structure 
and consequently also Tc's. 

Results: By numerical iteration of the mean particle 
number in the five angular momentum components and 
the s.c. mean-fields at fixed chemical potential we arrive 
at self-consistent solutions. For all plots the system size 
is 100 X 100 with at least ten sampling points per unit 
shift in particle number. Figure (3) shows the energy 
gap 2A, the norm of the s.c order parameter (defined 
as VTrOOt with O = El,s^^'^) at T=0 and as 
a ftinction of doped holes (up to 7 holes) for parameters 
u = -09, g = .06, A = .02 eV, both with and without 
Hartree-Fock terms. We find that Tg scales roughly lin- 
early with A and the reduced gap ^ ~ 3.2 is close to 
the weak coupling BCS value (3.53). For the calculation 
with H-F terms the magnitude of the order parameter 
(not in the figure) fits very well with 2 A. Without H-F 
terms there is a deviation from this fit around 5 holes, 
due to a momentum dependent gap. 

Fig. 4 shows the density of states, for the same pa- 
rameters and with Hartree-Fock terms and the DOS from 
the noninteracting band structure (8), i.e. the DOS with- 
out H-F terms. Since we find self-consistent solutions of 
both the band fillings and the gaps we calculate the DOS 
at finite temperature, above , for the realization with 
H-F terms. 



FIG. 4. Density of states. The thick solid line is for pa- 
rameters as in fig. (3), calculated at 150K. The thin (dashed) 
line is the DOS from band structure at T=0 (T=150K). 

The values of A and u chosen here correspond roughly 
to the pcrturbative spectrum (fig. 2) sd U = 5, but 
we have reduced g significantly in order for the pairing 
attraction not to dominate completely. One could ar- 
gue that the second order perturbation theory can be ex- 
pected to overestimate the core polarization effect that 
gives the 5^ term. Of course, the actual magnitude of 
Tc and the gap that we find should not be taken too lit- 
erally since parameters are only rough estimates and we 
are doing mean-field theory at relatively strong coupling. 
The important result is the qualitative behavior of Tc as 
a function of doping. 

Without Hartree-Fock terms, Tc follows the DOS of 
the band structure with a corresponding maximum at 
the half-filled band. This should be contrasted with the 
results for the full Hamiltonian which has a Tc that is 
maximized close to three holes and which is strongly sup- 
pressed at five holes. The DOS is suppressed at four 
and six holes due to the commensurate lock-in discussed 
above. In fact, the influence of these special fillings is 
such that the DOS is low in the whole region between 
four and six holes where it is the highest without inter- 
actions. 
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FIG. 3. Spectral gap 2A and norm of Osc at T=0 and Tc 
for parameters u = .09, g = .06, A = .02 eV as a function 
of doped holes. The lower (upper) curves are with (without) 
Hartree-Fock terms. 
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